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ABSTRACT 

A detailed record of the physical processes that operate during post-main-sequence evolution 
is contained in the internal chemical structure of white dwarfs. Global pulsations allow us to 
probe the stellar interior through asteroseismology, revealing the signatures of prior nuclear 
burning, mixing, and diffusion in these stars. I review the rapid evolution of structural models 
for helium-atmosphere variable (DBV) white dwarfs over the past five years, and I present 
a new series of model-fits using recent observations to illustrate the relative importance of 
various interior structures. By incorporating physically motivated C/O profiles into double- 
layered envelope models for the first time, I finally identify an optimal asteroseismic model 
that agrees with both diffusion theory and the expected nuclear burning history of the pro- 
genitor I discuss the implications of this fundamental result, and I evaluate the prospects for 
continued progress in the future. 
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1 INTRODUCTION 

The internal structure of a white dwarf star contains a detailed 
record of the physical processes that operate during post-main- 
sequence evolution. The composition of the deep interior is de- 
termined by the nuclear reaction rates during core helium burning 
in the red giant progenitor, while the chemical profile in the outer 
core is set by helium shell burning and various mixing processes 
including convective overshoot. More than 99% of the mass in a 
typical white dwarf is contained in the C/O core, which is embed- 
ded in an envelope of C, He, and H with traces of heavier elements. 
Near the end of post-AGB evolution the white dwarf envelope can 
become uniformly mixed during a very late thermal pulse, which 
bum s off most of the residual H and creates a bom-again AGB 
star ilben et aljlT983l) . As this newly formed H-deficient star be- 
gins to descend the white dwarf cooling track, the high surface 
gravity (log ^ ~ 8) coupled with chemical and thermal diffusion 
leads to compositional stratification. The lighter elements float to 
the surface, slowly building a nearly pure surface He layer above 
the remainder of the envelope as the star cools. 

Asteroseismology is the only observational method avail- 
able to probe the interior structure of white dwarfs and cali- 
brate theories of nuclear burning, mixing, and diffusion in these 
stars. The helium-atmosphere (DB) white dwarfs are particularly 
suitable for asteroseismic study, since they are observed to pul- 
sate in a range of effective temperatures between 22,400 and 
28,400 K (jBeauchamp et al. 1999). The mass of the pure He sur- 
face layer is predicted to grow by roughly an order of magni- 
tude within this range, so an asteroseismic calibration of diffu- 
sion theory is possible by studying DB variables (DBVs) at var- 
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ious temperatures ('Dehner Sl Kawaler" 199^; 'Fontaine & BrassarJ 
I2OO2; Althaus & Corsico 2004 ; Metcalfe et al. 2005) . Unlike the 
hotter class of variable white dwarfs (the PG 1159 stars) DBVs do 
not experience significant gravitational contraction, which simpli- 
fies the asteroseismic anal ysis since the m echanical and thermal 
structure are decoupled 1 Winget et al J2 00'7). The cooler hydrogen- 
atmosphere variables (DAVs) also enjoy this advantage, but their 
cores are significantly more degenerate — reducing the influence of 
the deep interior structure on their pulsation periods. 

Despite their extraordinary potential, only two DBV white 
dwarfs have been sufficiently characterized to allow detailed astero- 
seismic study. The brightest member of the class, GD 358, has been 
the target of three multi-site observin g campaigns of the Whole 
Earth Tele scope (WET; Nather et al. 1990; "Wing et et alJ Il994l; 
IVuille et' al. 2000; Kepler et al. 2003). These data revealed a series 
of 1 1 dipole modes (£=1, )!7=0; Kotak et al. 2003) of consecutive ra- 
dial overtone (fc=8-18; .Bradlev & Wingeliil 994) that could be used 
for model-fitting. Prior to the Sloan Digital Sky Survey (SDSS; 
iNittaet ali i2005) the faintest member of th e class was CBS 1 14 , 
which was first subjected to detailed study b vlHandler et alj j2002h 
A recent dual-site campaign on this star bv lMetcalf^^] <2005h 
documented a total of 1 1 dipole modes with radial overtones rang- 
ing from k=^ to k=2Q. 

In the following section, I review the recent development of 
detailed structural models for DBV white dwarfs with an empha- 
sis on the published models for GD 358 over the past five years. 
I then present a series of new model-fits to the latest observations 
of CBS 1 14 in Section|3|to illustrate the relative importance of the 
various interior structures expected from stellar evolution theory. 
The final fit in this series represents the first successful application 
of a model containing a complete, physically motivated description 
of the white dwarf interior. Finally in Section|4| I discuss the impli- 
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cations of this result for theories of diffusion and nuclear burning, 
and I evaluate the prospects for future progress. 



2 MODEL DEVELOPMENT 

Although detailed observations of the DBV white dwarf GD 358 
have been available since 1994, the realization of their full po- 
tential had to wait for the large-scale exploration of models made 
possible by fast, inexpen sive computers. The original analysis by 
iBradlev & Winged il994) attempted to match the observed pulsa- 
tion periods using a small grid of models around an initial guess 
based o n general scaling arguments and analytical relations devel- 
oped b vlKawalei* h99ab.lKawaler & Weiss! Jl990l) . lBrassard et alj 
il992h . and Bradle \^^iT ' ^^3l) . The main physical parameters 
that were adjusted to fit the models to the data included the stellar 
mass (M»), the effective temperature (Ten), the mass of the surface 
He layer (Mm), and several fixed core profiles including uniform 
and variable C/0 mixtures. The best-fit model from this analysis 
matched the observed periods with a precision of ap = 2.69 s, but 
had a surface He layer (Mhc ~ 10"^) that wa s thinner than expected 
from stellar evolution theory (Mhc ~ 10~^: IPantona & Mazzitellil 
Il979h . 

This pioneering work guided the development of a 
new model-fitting method based on a parallel genetic al- 
gorithm ( Metcalfe & CharbonneaJ 1200 3) running o n a 64- 
node commodity-hardware Linux cluster (Metealf^^Nathed 
| 200Cl). Using essenti ally the same models as iBradlev^^^ngeg 
iMetcalfe et alJ J2000h had the genetic algorithm explore a much 
broader range of the three main parameters (M, , rcff,MHc) with 
a similar set of fixed core profiles. This innovative search re- 
vealed a second family of solutions with the expected thick surface 
He layers, outside of the range considered by the original study, 
which provided a considerably better match to the observations 
(ap = 1.50 s). The results also demonstrated that the models were 
much more sensitive to the core composition than was previously 
believed. The conventional wisdom held that ^-mode pulsations in 
white dwarf stars were primarily envelope modes. While it is true 
that the horizontal displacements in the envelope from the non- 
radial oscillations are so large that they make the motions in the 
core almost invisible, the inner 90% of the mass actually contains 
the first several nodes of the radial eigenfunction for the pulsations 
observed in DBV stars. 

Motivated by this newly-discovered sensitivity to the core 
composition, M etcalfe et al. (2001) extended the genetic algorithm 
fitting method to optimize two additional parameters describing 
a generic C/0 profile. The adjustable central oxygen mass frac- 
tion (Xo) was fixed to a constant value out to some fractional 
mass (q) where it then decreased linearly in mass to zero oxy- 
gen at 0.95 Mr/Mf The optimization of these five parameters 
led to a model that matched the observations with a precision of 
(Tp = 1.28 s, representing a substantial improvement over the best 
strictly 3-parameter model (ap = 2.30 s). From this first asteroseis- 
mic measurement of Xn. iMetcalfe et all |2001) attempted to con- 
strain the rate of the astrophysically important '^C(Qf,7)'*0 nu- 
cl ear reaction, extra polating from the publishe d chemical profiles 
of'Salarisetal.'('l997). A follow-up study by Metcal fe et all <2002h 
quantified many possible sources of systematic error in this mea- 
surement, but failed to identify any that were substantially larger 
than the statistical errors. However, persistent disagreement be- 
tween the derived values of Teff and M, and the spectroscopic esti- 
mates suggested room for improvement in the model. 



Ipontaine & Bra ss^S i2002h disputed the results of 
Metcalfe et al. (2001), arguing that the model envelopes should 
include a double-layered structure. If we assume that DBV white 
dwarfs are the evolutionary descendants of the hotter PG 1159 
stars, then the surface He layer should be only Mue ~ 10"* M, , 
situated above the still-uniform envelope (Mem ~ 10"^ M*). This 
could explain the tw o families of solutions for Mhs found by 
Metcalfe et al. ( 2000) when assuming a single-layered structure. 
To illustrate the po tential of double-layered envelope models, 
iFontaine & BrassarS attempted to match the pulsation periods 
of GD 358 with a targeted grid. Using several uniform core 
compositions, they adjusted four parameters (M, , TcffjMcnv.MHc) 
and found a best-fit model with ap = 1.30 s. This suggested that 
double-layered envelope models with no structure in the core 
could explain the observat ions of GD 35 8 nearl y as well as the 
single-layered models of IMetcalfe et alj i200ll) which had an 
adjustable C/0 profile in the core. 

Around the same time, iMetcalfel i2003h identified two large 
sources of systematic error in the determination of Xq. Previous 
fitting had attempted to match not only the periods, but also the 
forward period spacing (AP = P^+i - Pi,) simultaneously. Ironi- 
cally, this was intended to minimize the impact of systematic er- 
rors on the calculated periods, but it biased the determinations 
of Xo to be higher. A second bias in the same direction came 
from fixing the mixing-length/pressure scale height ratio too high, 
leading to overly efficient convection that modified the thermal 
structure of the models. By matching only the periods and usirig 
the ML2/a=1.25 prescript ion for convection I Bohm & Cassinelli| 
1971; Bea uchamp et a l. 1999), the 5-parameter single-layered en- 
velope fit for GD 358 improved to ap = 1.05 s and the derived value 
of Xo suggested a '^ C(a,7)'*0 rate that was consistent with labo- 
ratory measurements jAngulo et aljl999t) . 

It finally became clear that neither of these models repre- 
sented an adequate description of the white dwarf structure when 
iMontgomerv et all i2003h identified an inherent symmetry in the 
way the pulsations sample the interior of DBV models. They found 
that generic perturbations caused by structure at specific locations 
in the core could not be distinguished from structure at correspond- 
ing locations in the envelope. This core-envelope symmetry could 
produce confusion between a C/0 gradient near 0.5 from 
prior nuclear burning, and the base of the surface He layer in the 
outer 10"* M» from diffusion. However, by using realistic internal 
chemical profiles in both regions of the models, they determined 
the relative importance of the structures expected at various loca- 
tions. This led them to conclude that the deep interior core structure 
should leave the largest imprint on the observed periods, while the 
base of the surface He layer was expected to have the smallest in- 
fluence. The total envelope mass was less important than the core 
structure, but had a larger effect than the outermost layer. 

To investigate whether the str ucture in the core and envelope 
could be measured simultaneously. 'Metcalf e et alJi2003h modified 
their models to include a parameterization of the double-layered 
envelopes used by IPontaine & Brassa rd (2002). They perform ed 
4-parameter fits comparable to those of Fontaine & Brassard as 
well as 6-parameter fits that also included an adjustable C/O pro- 
file. Although the addition of core structure led to a significant im- 
provement, neither fit reached the level of precision attained by us- 
ing single-layered envelope models, and the disagreement with the 
spectroscopic estimates of Tcff and M, became even wors e. How- 
ever, by fitting the peri ods of the model from Fontaine & Brassar j 
i2002h . lMetcalfe et al J discovered a systematic offset in the temper- 
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Table 1. New model-fitting results for CBS 1 14. 





T,if (K) 


M*/M0 


log(Me„v/M.) 


log(MHe/M.) 


Xo 


1 


crp (s) 


3 


25,300 


0.675 


-2.80 








2.44 


4 


25,800 


0.630 


-2.42 


-5.96 






2.33 


5 


23,800 


0.690 


-2.85 




0.93 


0.47 


1.81 


6« 


25,200 


0.625 


-2.40 


-5.94 


0.91 


0.42 


1.51 


6b 


24,900 


0.640 


-2.48 


-5.94 


0.71 


0.38 


1.27 



" Generic C/0 profile from lMetcalfe et ai]l200ll) 

Physically motivated C/0 profile based on ISalarisetalJ<1997h 



ature scales that could be at tributed to the different radi ative opaci- 
ties used in the two models I'F ontaine & Brassar(il994) . 

Confronted with so many fitting results for one set of 
observations — some of them from models with very different inte- 
rior structures — a casual observer might lose all confidence in our 
ability to do asteroseismology. Ultimately this is a symptom of the 
fact that we are confined to forward modeling: we adopt a structural 
model and see how closely it can match the observations. We can 
inform this choice with theoretical expectations, but we will always 
find an optimal model — even when the model is incomplete. The 
results from an incomplete model are not meaningless; they provide 
clues that help us identify the missing ingredients. The two fami- 
lies of solutions for the surface He layer in single-layered models 
told us that the real star showed evidence of composition transition 
zones at two locations in the envelope. The systematic differences 
between the derived masses and effective temperatures compared 
to spectroscopic estimates hinted that we needed to update the con- 
stitutive physics. How closely we can match the observations with 
each of the incomplete models gives us insight about which struc- 
tural ingredients are the most important. 



3 NEW MODEL-FITTING RESULTS 

Recent observations of CBS 114 bv lMetcalfe et alj <2005h finally 
give us an opportunity to evaluate the models reliably on a second 
DBV white dwarf. If we fit these new data using each of the mod- 
els that have previously been applied to GD 358 — but keeping the 
model physics fixed — we can isolate the impact of various interior 
structures on the quality of the match to the observations. This ex- 
ercise may be easier for CBS 1 14 since the intrinsic core-envelope 
symmetry will affect slightly different regions of the model, so we 
might have less difficulty distinguishing the signatures of the sur- 
face He layer and the deep C/O profile. It should also give us a 
clearer understanding of what is actually being measured by each 
of the structurally incomplete models, and how the other parame- 
ters are affected by the incompleteness. 

All of the modeling discussed b elow was performed w ith the 
WD-40 code originally described bv lMetcalfe etalj l20(Xtl . with 
several important updates to the constitutive physics. The equa- 
tion of state (EOS) for t he c ores and the envelop es come from 
iLamb & Van Hornl il975h and lFontaine et alj il977h . respectively. 
To describe the composition trans ition zones in the envelope, we 
use an adaptation of the method o f'Arcor agi & Fontaine! iT98Cl) for 
single-layered models, or a parameterization of the time-dependent 
diffusion results of Fontaine & Brassard ( 2002) for double-layered 
models, as described bv .Metcalfe et aL li2003.) . The C/0 profiles in 
the core are set either by the generic parameterization developed by 
iMetcalfe et alj i200lh . or by a new version that uses an adjustable 



profile based on the calculations of lSalaris et al] jl 997|). The mod - 
els use the updated OPAL opacity tables ( lelesias & Rogers" 199^, 
neutrino rat es from I toh et al. ( 19962, ^^'^ mixing-length pre- 
scription of lB6hm'& C assinelli ( jlgTlj) with ML2/a=1 .25 convec- 
tive efficiency as suggested by B eauchamp et alj il999l) . 

The new set of optimal models for CBS 1 14 are listed in Ta- 
bleQ In each case, a model-fitting method based on a parallel ge- 
netic algorithm was used to search a broad range for each parameter 
and o ptimize the match between the observed and calculated peri- 
ods ( Metcalfe & Charbonneau 2003). The search range for each pa- 
rameter and the details of the model-fitting procedure have b een de- 
scribed by Metcalfe et al. ( 2000) for Np = 3.lMetcalfe et al.| <200ll) 
for Np = 5, and Metcalfe et aU t2003l) for Np=4 and A',, = 6. 

Beyond the basic stellar parameters (7^ff,M,), the simplest of 
the models (Np = 3) includes an adjustable mass for a single sur- 
face layer situated above a pure C core. The derived value of Teff 
falls comfortably within the range determi ned from spectroscop y 
(Teff = 23, 300-26,200 K, log^ = 7.98-8.00; Beauchamp et alJl999^ ■ 
but the derived mass (logg ~ 8.12) is significantly above the spec- 
troscopic estimate. Despite the many simplifications inherent in 
this model, the derived envelope mass is in the range expected 
from stellar evolution theory (Dantona & Mazzitelli 1979). As with 
GD 358, a second family of models with envelope masses near 
10"* M, produce better than average fits to the observations of 
CBS 114, suggesting that double-layered models might be more 
appropriate (see Fig.0. 

The optimal parameters for the double-layered envelope 
model (Np = 4) include values of Tcff and M, that are both consis- 
tent with the spectroscopic estimates, but the slight improvement 
to the fit is not sig nificant. According to th e Bayes Information 
Criterion (BIC; see^Mont gomerv et alJl200lh . the addition of one 
adjustable parameter to the fit should reduce the residuals from 
o-p = 2.44 s for A^p = 3 to crp = 2.19 s for Np = 4. The actual fit 
yields crp = 2.33 s for Np = 4, which is statistically worse than the 
single-layered fit. In the absence of core structure, single-layered 
envelope models provide a statistically better match to the obser- 
vations of CBS 114 than double-layered envelope models. Similar 
results were obtained for GD 358, where the double-layered fit of 
IMetcalfe et Z|j2003l ctp = 2.17 s) was also statistically worse than 
the single-layered fit of Metcalfe et al. ( 2000, crp = 2.30 s) for pure 
C cores. Even so, the derived values of Mc„v and Mue from both 
stars qualitatively ag ree with the expectations of diffusion theory 
jMetcalfe et alj200a) . suggesting that the double-layered fits might 
be measuring real structure even if the models are incomplete. 

By contrast, the addition of an adjustable C/0 profile to single- 
layered models (A',, = 5) leads to a dramatic improvement in the 
residuals. The BIC leads us to expect a reduction from ap = 2.44 s 
to crp = 1.96 s, just from the addition of two adjustable parame- 
ters. In fact, the single-layered fit with a C/0 core improves the 
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Figure 1. Front and side views of tlie search range for: a the single-layered 
envelope [Np = 3) fit, and b the double-layered envelope (Np — 4) fit to 
CBS 114. Models that produce residuals within Icr (red), 3(t (yellow), lOcr 
(green), 25(t (cyan), and 40ct (blue) of the optimum are shown. 

fit to crp = 1.81 s, a statistically significant reduction. However, the 
derived values of Teff and M, do not agree nearly as well with spec- 
troscopic estimates, and the value of Xo is difficult to reconcile with 
the expectations from prior nuclear burning. Again , similar results 
were found for GD 358 bv iMetcalfe et al.l i200ll) . Despite these 
systematic errors, the derived values for Mcnv from the two single- 
layered fits agree with each other, suggesting that thick envelopes 
are a qualitatively robust structural feature of the results. 

For models that contain an adjustable C/0 profile in the core, 
the inclusion of double-layered envelopes (Np = 6) produces a sig- 
nificant improvement in the residuals, and restores the agreement 
with spectroscopic estimates of T^ff and M, . The residuals from the 
A'p = 5 fit should decrease from crp = 1.81 s to ap = 1.62 s from 
the addition of one adjustable parameter, but the A'p = 6 fit with a 
generic C/O profile actually decreases the residuals to crp = 1.51 s. 
This statistically significant improvement is accompanied by de- 
rived values of M^m and Mh^ that once again qualitatively agree 
with diffusion theory. The only remaining concern with this model 
is the implausibly large derived value of Xq. 

This may not be too surprising, considering the simplistic 
form of the generic C/O profile compared to chemical profiles gen- 
erated by realistic stellar evolution calculations. Although a con- 
stant value of Xo out to som e fractional mass (q) is a generic fea- 
ture of such simulations (see'Metcalfe et all2003l their Fig. 2), the 
O profile in the outer core is expected to be m ore complicated than 
a linear decrease. For example, the models of ISalaris et alj il997l) 
show an initially sharp drop in the O mass fraction near the edge 
of the progenitor's convective core, surrounded by a more gradual 
decline produced during He shell burning, followed by an abrupt 
drop where the temperature and pressure in the progenitor were no 
longer sufficient to sustain shell burning. A new implementation 
of the adjustable C/O profiles continues to keep Xo constant out to 
some fractional mass q, but then falls to zero oxygen at 0.95 M,-/Mf 
with a more physically motivated shape. ForXo = 0.76 and q = 0.46, 




0.2 0.4 0.6 0.8 



M/M. 

Figur e 2. A comparison of the generic C/O profile from |Metc^f^^^ 
(2001, red) with the physically motivated profile based on Tsal^^^L 
( 1997, blue) used for the two N,, = 6 fits to CBS 114. The plotted values 
of Xq and q represent the final entry in TablefTI 

the profile closely resembles the results of lSalaris' etalJjl99l their 
Fig. 3), and for any other values of Xo and q the profile is scaled to 
maintain a similar shape between q and 0.95 Mr/M, (see Fig.|2j- 

By incorporating these physically motivated C/O profiles 
into double-layered envelope models, the optimal parameters for 
CBS 114 lead to an interior structure that agrees with both diffusion 
theory and the expected nuclear burning history. Without adding 
any adjustable parameters, this modification to the C/O profile in 
the outer core reduces the residuals from crp = 1.51 s to crp = 1.27 s, 
an improvement comparable to what would be expected from 1.5 
additional free parameters. The derived values of Tcff and M» are in 
reasonable agreement with spectroscopic estimates, the structure of 
the envelope is approximately as expected from diffusion theory, 
and the core C/O profile implies a rate for the '^C(a,7)"'0 reac- 
tion that agrees with recent laboratory measurements. We should 
ultimately be able to refine theories of diffusion and nuclear burn- 
ing to accommodate these new observational constraints. 



4 DISCUSSION 

The past five years have seen unprecedented progress in the devel- 
opment of white dwarf interior structure models, and in our ability 
to fit them to the available observations. This rapid improvement 
in our understanding has been most evident for the DBV white 
dwarfs, where the physical conditions are ripe for asteroseismic in- 
vestigation. Our models for these stars have evolved from a simple 
parameterization of a single surface He layer above a pure C core 
to the double-layered envelope structure expected from diffusion 
calculations surrounding a physically motivated C/O profile in the 
core. We have finally been able to demonstrate that these structural 
ingredients, when taken together, lead to significant improvements 
in our ability to match the observations of CBS 1 14. 

It is now clear that the detailed core C/O profile is the most 
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important feature for quantitative asteroseismology of DBV white 
dwarfs. The most dramatic improvements in our ability to match 
the observations came from the addition of C/0 cores, even when 
using relatively crude approximations for th e structure in the enve- 
lope. The generic C/0 profiles developed bv lMetcalfeetalJj200ll) 
seemed to reproduce the most important asteroseismic features, but 
the physically motivated profiles introduced here appear to be nec- 
essary for a reliable determination of the composition deep in the 
core. Disentangling the roles of the '^C(a,7)'^0 reaction and mix- 
ing processes such as con vective overshoot shou ld be possible us- 
ing the derived value of q iStraniero et allZOOBl) . New simulations 
of white dwarf interior chemical profiles that attempt to match the 
asteroseismic C/0 profile will help to calibrate these processes. 

Without a C/0 profile, the double-layered envelope structure 
expected from diffusion actually diminishes the ability of the mod- 
els to match the observations, compared to simpler envelopes. We 
can now understand this to be a natural outcome of the relative 
importance of various regions of the interior: C/0 profiles leave the 
largest imprint on the pulsation periods, followed by the base of the 
envelope (Mcnv) and the surface layer (Mne). In the absence of core 
structure, models that maximize the imprint of Menv provide better 
fits. Since double-layered models effectively split the He transition 
zone between two separate locations in the envelope, the imprint 
of Mcnv is reduced relative to single-layered models. With the C/O 
core included, the fine tuning made possible by the double-layered 
models produces a significant improvement and leads to a credible 
envelope structure. Future time-dependent diffusion calculations 
should attempt to reproduce these asteroseismic measurements. 

We can instill greater confidence is these results by apply- 
ing our physically motivated structural models to additional white 
dwarf stars. Imposing the requirement that our model must be able 
to explain two independent sets of observations — from stars oper- 
ating under slightly different physical conditions — is a very pow- 
erful constraint. Unfortunately, GD 358 may not be the best choice 
for such an experiment because the intrinsic core-envelope sym- 
metry in this case appears to create an intractable ambiguity be- 
tween the locations of the surface He layer and the core C/0 pro- 
file, even when using physically motivated shapes for the chem- 
ical transitions. A better candidate might be found among DBV 
stars that are hotter than CBS 1 14 (e.g. EC 20058-5234) or cooler 
than GD 358 (e.g. PG 1456+103). Many new DBV st MS have re- 
cently been identified in the Sloan Digital Sky Survey il Nitta et all 
120051) ■ but most of them are considerably fainter than the previously 
known sample. Extensive follow-up observations to identify the 
most promising multi-mode pulsators will require 2-m class tele- 
scopes e quipped with time-series o ptimized frame-transfer CCD 
cameras jNather & Mukadamll2004h . With the computational tools 
in place, we should quickly be able to refine our understanding of 
the forces that shape white dwarf interiors. 
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